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Charged particle pseudorapidity distributions are presented from the PHOBOS experiment 
at RHIC, measured in Au+Au and Cu+Cu collisions at ^/s NN =19.6, 22.4, 62.4, 130 and 
200 GeV, as a function of collision centrality. The presentation includes the recently analyzed 
Cu+Cu data at 22.4 GeV. The measurements were made by the same detector setup over a 
broad range in pseudorapidity, < 5.4, allowing for a reliable systematic study of particle 
production as a function of energy, centrality and system size. Comparing Cu+Cu and 
Au+Au results, we find that the total number of produced charged particles and the overall 
shape (height and width) of the pseudorapidity distributions are determined by the number 
of nucleon participants, N var t- Detailed comparisons reveal that the matching of the shape 
of the Cu+Cu and Au+Au pseudorapidity distributions over the full range of n is better for 
the same N P art/2A value than for the same N par t value, where A denotes the mass number. 
In other words, it is the geometry of the nuclear overlap zone, rather than just the number 
of nucleon participants that drives the detailed shape of the pseudorapidity distribution and 
its centrality dependence. 



I. INTRODUCTION 



The charged particle pseudorapidity distribution, dN c h/ drj 1 is a well defined experimental quan- 
tity that reflects the initial conditions of the system, e.g. parton shadowing and gluon saturation, 
and also the effects of rescattering and hadronic final state interactions: it represents the time- 
integral of the particle production throughout the entire collision. With the advent of Cu+Cu 
collisions at RHIC, the system size dependence of important observables can be studied using 
different collision geometries. The Cu+Cu results [l[ test the simple scaling features observed pre- 
viously in Au+ Au collisions 0, Q • They significantly extend the N par t range measured in Au+Au 
collisions, while the two systems can also be compared at the same N part , as illustrated below. 

II. EXPERIMENTAL SETUP AND DATA ANALYSIS 

The Cu+Cu data were collected with the multiplicity array of the PHOBOS detector Q during 
the RHIC 2005 run. The array consists of single-layered silicon sensors assembled into a long, tube- 
shaped Octagon detector surrounding the collision point, and into three Ring sensors on each side, 
detecting large- 1 77 1 particles. Simulations of the detector performance were based on the HIJING 
event generator and GEANT, including the response of the scintillator Paddle trigger counters. 

Data from the Cu+Cu and Au+Au collisions were analyzed using the 'hit-counting' and 'analog' 
methods The latter was corrected for multiply-charged fragments emitted at large ry. This 
correction decreases with centrality and collision energy, and it is less than 6% of the total number 
of charged particles. 

The estimated trigger efficiency is 84+5% and 75+5% in Cu+Cu collisions at 200 and 62.4 GeV. 
The centrality of the collision was estimated from the Paddle scintillator signals. At 22.4 and 
19.6 GeV, the pathlength-corrected energy sum [3[ deposited in the Octagon was used (\r]\ < 3.2). 
A Glauber- model calculation was employed to estimate (N part ) for each centrality bin. 

III. RESULTS 

The dN c h/drj distributions in Cu+Cu collisions for various collision energies and centralities are 
shown in Fig. [TJ On the right panel the Cu+Cu and Au+Au collisions are compared, where the 
centrality bins are chosen such that (N part ) in both systems are similar. One can conclude that 
although the distributions agree at the same (N part ) to first order, there are differences at large 
\r]\ and low energies. In this context, we note that the two nuclear spectator remnants are larger 
in Au+Au than in Cu+Cu collisions and thus may be playing a role. 

The dNch/drj distributions exhibit longitudinal scaling when observed from the rest frame of 
one of the colliding nuclei. The coordinate transformation to the 'target' frame approximately 
corresponds to a shift by the beam rapidity, ybeam- Figure [2] compares the dN c h/dr)' distributions 
(where 7/ = 77 — ybeam) after normalization by N part : a) data from the Cu+Cu and Au+Au systems 
plotted at the same fraction of the total cross section (0-6% most central bin), and b) at the same 
value of N part /2 A (where A is the mass number). Both cases indicate that the scaled particle 
density only depends on the collision energy and geometry, but not on the size of the nuclei. The 
dN c h I drj' I (Np ar t) distributions for the same centrality in both systems agree within errors, and the 
overall agreement improves if the centralities are compared on the basis of the N part /2A quantity 
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FIG. 1: Left panel: pseudorapidity distributions of primary charged particles from Cu+Cu collisions 
at 22.4, 62.4 and 200 GeV collision energy per nucleon pair for various centrality bins. Right panel: 
comparison of dN c h/dr) distributions in Cu+Cu and Au+Au collisions - corresponding to the same {N par t)- 
90% C.L. systematic errors are shown as bands. 



(the fraction of participating nucleons). The longitudinal scaling is similarly present in the Cu+Cu 
and in the Au+Au data. 

The factorization between collision energy and centrality can be most precisely studied by exam- 
ining the ratios of the dN / drj' / (N part ) distributions in central and semi-central collisions, denoted 

by RpQ rt , at various energies. The published Au+Au results @ are shown by the inset of Fig. [3^, 
exhibiting the same factorization feature as the recent Cu+Cu data. The above ratio for Cu+Cu 
data is similar to that in Au+Au data, except at the highest rf values. Fig. [5Jd shows the Rpc 
ratio for Cu+Cu and Au+Au for centrality bins where the N part /2A values are matched. The lat- 
ter quantity characterizes the initial geometry more precisely, and indeed, the centrality evolution 
of the dNch/dt] distributions measured in Cu+Cu and Au+Au collisions are most similar if the 
centrality is quantified by N part /2A. 

The total number of charged particles, N to t, normalized by N part , is presented in Fig. 0] for 
Cu+Cu collisions at 22.4, 62.4 and 200 GeV collision energy per nucleon pair, and compared to 
smaller (p+p, d+Au) and larger (Au+Au) systems as a function of centrality. One can conclude 
that Ntot scales approximately linearly with N part , and the normalized yield has similar values for 
the two heavy colliding systems. The d+Au data do not seem to interpolate smoothly between 
the p+p and heavy ion data points. 
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FIG. 2: Pseudorapidity distributions in Cu+Cu and Au+Au collisions at various RHIC energies, normal- 
ized by the number of participant pairs, plotted as a function of r{ = r\ — ybeam, for a) the 6% most central 
events and b) central events with similar N part /2A. 90% C.L. systematic errors are shown for a few typical 
data points. 




FIG. 3: The semi-peripheral to central dNch/drf / (N part /2) ratio for Cu+Cu and Au+Au collisions at 
RHIC energies. The centrality bins are selected a) according to fractional cross section (35-40%/0-6%) 
and b) such that Npart/2 A is matched between the two systems. Inset: the same quantity plotted for 
Au+Au data for four different energies. The error bars represent 90% C.L. systematic errors on the ratio. 



IV. SUMMARY 



Charged particle r\ distributions were presented, including the recently analysed Cu+Cu data 
taken at various collision energies. The dN/dr)' distributions scaled by N part /2, as well as their 
peripheral to central ratio were found to be independent of the mass number, A, of the colliding 
nuclei if centrality classes with the same N par t/2A (fraction of participant nucleons) are compared. 
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4: The integrated number of charged particles, scaled by N part /2, in p+p, d+Au, Cu+Cu and 



3d pa 

Au+Au collisions as a function of centrality 0,13 
bars. 



The uncertainty of N par t has been included in the error 
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